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Continuum Thermodynamics of the GluoA^c Plasma 

Saumen Datt£0 and Sourendu 

Department of Theoretical Physics, 
Tata Institute of Fundamental Research, 
Homi Bhabha Road, Mumbai 400005, India. 

We study the thermodynamics of SU(A'^c) pure gauge theories for A'^c = 3, 4 and 6. The continuum 
and thermodynamic limits of bulk quantities such as the pressure (p), energy density (e) and the 
entropy density (s) are taken by using several different lattice spacings and volumes. There is no 
window of temperature in which a non-trivial conformal theory describes bulk thermodynamics. We 
extract the latent heat of the first-order deconfinement phase transitions and observe good scaling 
with Nc- For all quantities that we measure, strong Nc scaling holds, except, possibly, very close to 
the transition temperature, T^; however we are unable to find strong evidence for scaling with the 
't Hooft coupling in thermal quantities at the small values of A'^c which we study. 

PACS numbers: 



I. INTRODUCTION 



Interest in SU(A^c) gauge theories with large Nc began with the pioneering studies of where it was shown 
that in 1 + 1 dimensions the limit Nc oo when taken such that the gauge coupling, g — > 0, with the 't Hooft 
coupling A = Ncg^ fixed gave rise to an interesting and non-trivial but tractable theory. In this so-called 't Hooft 
I scaling limit stable mesons exist with properties which parallel much of known hadron phenomenology. The limit is 
O-i non-perturbative, and the computation of any correlation function requires the summation of an infinite number of 
^ ' well-characterized Feynman diagrams. Corrections to this limit appear in powers of 1/Nc in general, and in powers 
of 1/N^ for the pure gauge theory. Since then many models have used this 't Hooft limit [2], including the currently 
fashionable conformal cousins of QCD: the so-called A/" = 4 supersymmetric theories which turn out to be tractable 
using the AdS/CFT correspondence. Lattice simulations of theories with small iVc > 3 can test the smoothness of 
, approach to the 't Hooft limit. 

■ A contemporary reason for studying the large-TVc theory is in the light it could throw on the phase diagram of 
Qv^ ] QCD. At any fixed Nc with two flavours of massless quarks, the theory is expected to have a second-order chiral 

■ symmetry restoring transition at finite temperature, T, which is in the 0(4) universality class. If a tiny mass is 
' given to the quarks, then the theory has a cross-over at finite temperature instead of a phase transition. The limit of 

infinite quark mass corresponds to the pure gauge theory, which has a first-order deconfining thermal transition for 
all Nc > 2. Hence at an intermediate quark mass there is a critical point, in the Ising universality class, which ends a 
T-H line of first-order deconfining transitions [sj. We represent this information in the diagram of Figure [TJ In the regions 
^ with the chiral cross-over (region A) and the large Nc deconfinement transition (region B) the phase diagrams in the 
, plane of T and baryon chemical potential, ^, are topologically distinct; a representative phase diagram from each of 

■ these regions is also shown. In region A one expects a first-order phase transition line dividing the chiral symmetry 
J-j [ broken phase from the deconfined gluoTVc plasma phase with a critical end-point at finite /x. In region B one might 

. - - 1 expect a triple-point with coexistence of the three phases — gluoiVc plasma, baryonic, i.e., chiral symmetry broken 
and confined, and quarkyonic, i.e., chirally symmetric and confined [J]. Since the arguments for the existence of a 
quarkyonic phase are based on a picture of the large- iVc theory which arises in the 't Hooft limit, it is interesting to 
explore the range of validity of such large- iVc arguments. 

In this paper our primary interest is in examining the thermodynamics of S\J{Nc) pure gauge theory in the large 
volume and continuum limits. The deconfinement transition in SU(Afc) theories has been studied in [^-Q- The latent 
heat of the transition was studied in '^,'8*1 on spatial volumes, V, of size < 3. The equation of state (EOS) has 

been studied with lattice spacin g of a — 1/(5T) up to 2Tc in [10] and up to 3.5Tc in [ll|. Previous work has shown that 
the continuum limit of the EOS |3l| is hard to extract on such coarse lattices In view of this, we have studied the 

equation of state at smaller lattice spacing in the extended temperature range O.QTTc T < ATc (preliminary results 
were presented in Q). We extrapolate to the continuum limit using the non-perturbative beta- functions determined 
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FIG. 1: The diagram on top shows the regions of first- and second-order thermal phase transitions for SU(A^c) gauge theories 
with two flavours of fermions as Nc and the quark mass, m, are varied. We also show the phase diagrams in two interesting 
regions — (A) the chiral cross over region and (B) the large- A^c deconfinement region. The topology of the phase diagrams in 
region A is expected to involve a line of first order transitions ending in a critical point at finite chemical potential ^. In region 
B large- A^c arguments indicate that the phase diagram contains three lines of first order phase transitions separating pairs of 
phases labelled baryonic, deconfined and quarkyonic, with a triple point where all three phases coexist. 

earlier [9|]. We perform finite size scaling studies by changing the spatial volume up to T\/V — 4 for a = l/(6Tc) and 
~ 3.5 for a = l/{8Tc). We use statistics significantly larger than used before in this context, and comparable 
to that used in studies of thermodynamics for SU(3) pure gauge theory. 

One of the purposes of a study like this is to perform lattice simulations for small Nc > 3, and, from measurements 
of any physical quantity, find the series expansion for it around Nc = oo in powers of 1 /Nc . In this work we assume 
that the Nc = oo limit exists and that there is a series expansion around it, and ask what our data imply for the 
radius of convergence of this series. Technically, this also means that for each series expansion, we ask how reliably 
the Nc oo limit can be taken from measurements for a few small values of Nc- 
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In agreement with earlier studies, we find good evidence for scaling to Nc oo keeping T/Tc fixed. Such a "strong 
Nc scaling", previously observed for many quantities on the lattice, also gives small l/N^ corrections. We also examine 
't Hooft's limiting procedure. It is clear that the notion of 't Hooft scaling has to be defined carefully in any theory 
with a non- vanishing beta- function ^2*1. Firstly because one has to use a (scale dependent) renormalized coupling 
with its attendant scheme ambiguities. Secondly because of this scheme dependence, l/N^ corrections may be moved 
between the operator expectation values and the 't Hooft coupling. Even keeping these uncertainties in mind, wc find 
that scaling at fixed A gives large corrections (in the expansion around Nc = oo) at small Nc, including the physically 
interesting case of TVc = 3. 

The plan of this paper is as follows. In the next section, we summarize the formula used for calculation of the 
various thermodynamic quantities. Next, in section Hill we discuss the latent heat of the deconfinement transitions. 
In section HVl we investigate the conformal symmetry breaking measure, which is the trace of the energy momentum 
tensor, A = e — 3p. In section |V] we discuss the pressure and the remaining bulk thermodynamic quantities. Section 
IVII is devoted to a comparison of results from theories with different number of colors, to get an estimate of the size 
of the leading 0{1/Nc) corrections. Section |VIl] analyzes the calculated thermodynamic quantities, to infer properties 
of the gluoTVc plasma. The appendix contains a detailed discussion of the beta-functions used in this study. 



II. FORMALISM AND DEFINITIONS 



The thermodynamics of the SU(A'c) gauge theory can be obtained from the partition function, 

Z{V,T)^ fvUcxpi-S), S = (5 J2 {l-RealP^,(2;)} (1) 

calculated on a space-time lattice with Ng lattice sites in each of the spatial directions and Nt in the time direction; 
the lattice sites are labelled by the 4-component index x and directions by Greek indices, 1 < /J., < 4. The bare 
gauge coupling, = 2Nc/ 13 determines the lattice spacing, a, which is implicit in the above equations. The spatial 
volume is = N^a^ and the temperature is T = l/{aNt). Since we investigate finite size effects, it is useful to 
introduce the aspect ratio, ( — Ng/Nt = T\/V . Pp,v[x) is the trace of the product of SU(A''c) valued link matrices 
\J around the plaquette in the /i, v plane starting at site x. The trace is normalized such that P^ti/(x) = 1 if the link 
matrices are set to the identity. 

The expectation value of the plaquette, 

E RealP,.(x), (2) 

over the ensemble at any temperature, T, is one of the primary observables on the lattice. The other is the expectation 
value of the Wilson line, 

1 ^* 

^ X x^ — l 

which is the spatial average of the product of time-like link variables wrapping around the lattice in the time direction 
at each spatial site, x. (L) is the order parameter of the confinement-deconfinement transition and changes from 
zero to non-zero values at the (first-order) phase transition temperature T^. In the deconfined phase (L) has 
allowed values. Often {\L\) is examined, although it is not an order parameter, since it has only two allowed values 
in the transition region: being close to vanishing in the confined phase and non-zero in the deconfined phase. These 
observables and the scaling of T^. to the continuum have been reported earlier Q. 

The integral in ([l]) is performed by Monte Carlo sampling, using a combination of pseudo-heatbath and over- 
relaxation steps, where all SU(2) subgroups of the SU(iVc) group elements are touched. For details of the algorithm 
used and its performance, see We studied SU(4) and SU(6) theories in the temperature range between Tc, and 
about 4Tc. Since a major focus of this study is to get results in the continuum and thermodynamic limit, over the 
whole temperature range, we use two different lattice spacings, a = 1/(6T) and 1/(8T), at each temperature and 
several different C The temperature scale for these theories was set in Q, where it was shown that near Tc 

the results are in the scaling regime. With the running of the coupling obtained in that study, we found very good 
agreement between the thermodynamic quantities extracted on lattices with the different a at all T. We have also 
made a few simulation runs for SU(3) gauge theories at large aspect ratios, to complement existing studies of latent 
heat in SU(3) gauge theories. 
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Bulk thermodynamic quantities are obtained by taking suitable derivatives of the partition function, Z(V,T). In 
particular, the energy density and pressure are given by 



T2 dlogZ 



V dT 



p = T 



dlogZ 



dV 



(4) 



The entropy density is given by the identity s/T^ — (e + p)/T^. 

A quantity that is easy to calculate on the lattice is the trace of the energy-momentum tensor, A = e — 3p. This 
is of some interest for models of the QCD plasma, since it is a direct measure of the conformal symmetry breaking. 
Using the above relations it is easy to show that 

^ = ^^t7^{SP), where 5P = P(/?, T) - (P(/3, T = 0)). (5) 
J * o log a 

The expectation value of 6P must be taken over the finite temperature ensemble. The subtraction of the plaquette 
expectation value at T = must be done at the same lattice spacing as the finite temperature simulation. This serves 
to remove an ultraviolet divergence from the plaquette. It also makes sure that A/T^ vanishes at T = 0, since both 
the pressure and the energy density vanish there. The derivative multiplying this non-perturbative factor is closely 
related to the beta-function of the theory. In the appendix we have a discussion of the choices of beta-functions and 
their influence on A/T^. 

The suggestion of [l^ was that since p = TlogZ/V, in the thermodynamic limit, the pressure can be evaluated 
with respect to some reference value by integrating the plaquette expectation value — 

P 1 1 7 Po , T f\^, dlogZ , ir.. 



Here po is the pressure at the reference temperature Tq = 1 / Nta{Po) . For the gluoiVc plasma the pressure is expected 
to be very small below T^, even as close to as T '-^ 0.9Tc. Conventionally one takes the reference temperature to 
be such a value and sets po = 0, as in the second expression, where the derivative has also been written out explicitly. 
Once A and p are known, e and s can be evaluated. 

Asymptotic freedom in SU(A^c) gauge theories lead us to expect that at sufficiently high temperatures one should 
obtain a free gas of gluons. Then thermodynamic quantities reach their ideal gas (i.e., Stefan-Boltzmann: SB) limits. 
There are lattice corrections to this limit 14]. When the pressure is evaluated through the integral method one has — 

^2 J . o 2 



PSB ^ d,A^,^^^ ^ _ ^^A7\ 1 , L 

^2 



^=3i-^== -G(Nt) where G(Nt) = l + t + (7) 



and (Ia — — 1. Since aT = ^/Nt, the correction terms come in powers of the lattice spacing a, and vanish in 
the continuum limit. The full G{Nt) is also known exactly from numerical computations, and listed in [l^ . When 
we discuss the numerical computations later we will use this full computation of G{Nt) and not the series expansion 
above. The difference between the two is about 2% for Nt = 6. We draw attention to the factor of (Ia- in the large 
Nc limit it is often replaced by N^, but at the small values of Nc we used, the difference is statistically significant. 
We use dA in this work, and thereby subsume some of the formally sub-leading 1/Nc corrections into this scaling. 

One of the pieces of physics we are interested in is the latent heat of the transition. In a thermodynamically large 
volume this is defined by the formula 

^ = lim _ e{T,-dT) \ ^ ^.^^ f AjT, + ST) _ AjT, - ST) \ 

where the second equality follows from the fact that p is continuous across a first-order phase transition. These 
formulae cannot be directly used at finite volume. We describe a method for determining Ae/T^ in section Hill 

Some remarks about the computation are best placed here. It was observed earlier that at lattice spacing of 
a < l/(6Tc) the two- loop beta- function with a l/N^ correction provides a good description of the change of measured 
length scales with the gauge coupling g^. Therefore, one should be able to perform continuum extrapolation of 
thermodynamic quantities for T > Tc using data on lattices with Nt > 6. This expectation should be correct except if 
there are large corrections in powers of a to the operators involved in defining components of the energy-momentum 
tensor. We show evidence later that there are no such large corrections. Another possible subtlety could reside in 
having to take the ^ — >■ oo (thermodynamic) limit before taking the continuum (a — >■ 0) limit. Such subtleties arise 
only when there are large correlation lengths, ^. Here we have a first order phase transition with ^ < 1/T [TBI. [l6|. 
Hence the full machinery of finite-size scaling need not be invoked if due caution is exercised: the extraction of the 
latent heat is one such case. We will mention these tests at appropriate places in the remainder of the paper. 
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III. THE LATENT HEAT 




FIG. 2: Histograms of \L\ and 5P normalized to have unit area under the curves. A muhi-peak structure in the histogram 
corresponds to muhiple (local) minima of the constrained free energy. As A'^s changes between the values shown, in steps of 2, 
the histograms move up or down monotonically in the region where they are labelled. 



When two copies of a system with a first-order transition are held at temperatures Tc ± 5T, the difference in their 
energy densities is finite. Since the transition is of first-order, this difference remains finite even when 5T — ^ 0, provided 
one has thermodynamically large systems. When V is large but not infinite, the correct limit may be obtained as long 
as 5T/Tc is much larger than 0{£,^ /V). If V is large enough, then this allows one to come close enough to Tc to have 
confidence in the result. For smaller S,'^ /V , this procedure breaks down. By examining the reason for the breakdown, 
we can develop a different procedure to estimate the latent heat. 

Finite-size effects at first-order phase transitions are best understood in terms of constrained free energies, F^. 
This is the free energy of a system in which \L\ or 5P are restricted to have fixed values but other parameters are 
allowed to vary according to the temperature. Fc near a first-order phase transition has multiple minima: the deepest 
corresponds to the value the constrained variable has in the stable phase when V ^ oo. As T changes, the depth 
of the minima change, and on crossing the deepest minimum flips. However, because the finite system sees finite 
barriers between the minima, the system explores all the phases. Hence the discontinuity is rounded off. As the 
thermodynamic limit, V — >■ oo, is approached, the "wrong phase" minimum becomes infinitely higher and the barrier 
separating it from the true vacuum also becomes infinitely high. As a result, the transition sharpens and gives the 
correct thermodynamic limit. 

The histogram of an observable obtained from its Monte Carlo history is proportional to exp{—VFc/T). The object 
of a finite-size scaling study of something like the latent heat is to be able to identify the thermodynamically stable 
phase from histograms such as those in Figure [H The same figure also illustrates the problem which is usually faced 
in computing the latent heat in SU(iVc) gauge theories. Although the histogram of \L\ has multiple well-identifiable 
maxima, the histogram of 6P has a single peak. If the specific heat in each of the phases is large, then the two- 
peak nature of A/T* could well be hidden until extremely large volumes are reached. Finite-size scaling methods 
were developed in the past to extract reliably the specific heat when multiple maxima are clearly developed (itI . 
However they are no t ap plicable here, and we need to use different techniques. We adapt one which was first applied 
to gauge theories in [l9j . 

Since the phases are well resolved by \L\, we can try to use the following criterion. The cold confined phase could 
be identified by requiring that \L\ < Lc, and the hot deconfined phase by \L\ > L^. The results of such a phase 
separation are stable as long as Lc and Lh both lie in the valley between the peaks of Examples of the probability 
density of SP obtained in the cold and hot phases so defined are shown in Figure [3] The figures illustrate the fact 
that these probability densities are very stable — at each fixed volume the pure phase probability distributions of SP 
are identical for a range of couplings around Tc- If we had not separated out the panels for different /3, the curves for 
the three different cases would have been indistinguishable apart from small statistical fluctuations. Furthermore, as 
one changes V (with fixed a) the mean value of 6P in the hot and cold phases remain the same. It appears that the 
difference between the values of 6P in the two phases is very stable under the variation of both T and V near Tc- It 
also turns out to be fairly stable under changes of Lc and L^. 

The agreement of the histograms of SP for the different /3 in the transition regime show that it is possible to reliably 
extract the limiting values of A/T"* in each of the phases. Examination of ([5]) shows that knowing the difference in 
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FIG. 3: Histograms of SP in the two phases of SU(4) gauge theory as defined using Lc = 0.024 and — 0.031. The panels on 
the left are results from 6 x 16^ lattices and those on the right from 6 x 20"^. The histograms are normahzed so that the area 
under each is unity. The pure phase histograms at three different couplings exactly overlap on each lattice. 



SP between the hot and cold phases one can extract easily the jump in A/T"* and hence the latent heat density, since 
the remaining factors are well-understood. In Table |T] we summarize our results for the latent heat for SU(Afc) gauge 
theories with Nc — 3, 4: and 6. In the entry for the latent heat, the first error is statistical while the second is a 
systematic error, i.e., the change in the result if — L^. changes by ± 20%. The results of are higher, but those 
of are consistent with ours, within the larger statistical and systematic errors of that study. 

In SU(4) gauge theory we found that the extracted value of Ae/T^ is fairly stable at fixed a as we change V. In fact, 
as is clear from Table HI when a — 1/ (6Tc) we found that there is no statistically significant change in the estimate 
of this quantity for ( > 3. The results for a = l/{8Tc) are consistent with this conclusion. For the SU(3) theory, on 
the other hand, it seems that C ^ 6 is needed for an estimate of the latent heat with equally small systematic errors, 
i.e., finite volume corrections are larger for the SU(3) theory. One sees that in going from A'c = 3 to 4 the latent heat 
density scales faster than cIa ■ This ties in with the intuition developed from a study of correlation lengths [H, [l^ 
that the SU(3) theory is weakly first-order. One expects that correlation lengths should also become shorter with 
increasing Nc iSi], and hence finite volume effects should be less pronounced. 

Following the analysis of the appendix, we understand that the lack of clear scaling of Ae/ {cIaT'^) to the continuum 
limit is not due to the use of an inappropriate beta-function. One possible explanation is that the large finite volume 
effects mask the approach to the continuum. If so, then one should be able to eliminate it by using another quantity 
with the same effect. Since A/T** has a very sharp peak as a function of T, and that measurement would also be 



7 



Nc Nt Ns P \L\c \L\h ^ejTt ^ej 1^^. 

16 5.6908 0.055 0.075 2.06(1) (3) 

4 24 5.6919 0.055 0.075 1.93(1)(3) 

g 32 5.6922 0.055 0.075 1.90(2) (2) 

24 5.8934 0.02 0.03 1.79(2)(4) 0.65(2) 

6 32 5.8938 0.02 0.03 1.54(2)(5) 

48 5.8940 0.022 0.032 1.44(4) (3) 

8 32 6.0609 0.013 0.019 1.67(4)(4) 0.68(3)' 

16 10.79 0.024 0.031 4.85(5)(6) 0.88(2) 

6 20 10.79 0.024 0.031 4.64(4) (5) 

4 24 10.79 0.024 0.031 4.57(4)(3) 0.85(2) 

22 11.08 0.013 0.017 4.58(5)(6) 

8 24 11.08 0.013 0.017 4.32(6)(6) 0.82(2) 

28 11.08 0.013 0.017 4.33(8)(6)) 

g 14 24.84 0.025 0.03 12.20(10) (4) 

6 18 24.84 0.025 0.03 12.47(4)(2) 0.92(2) 

8 20 25.46 0.012 0.015 11.93(34)(5) 0.90(3)' 



TABLE I: The latent heat of SU(A'^c) gauge theories for A^'c = 3, 4 and 6. The thermodynamic limit of Ae/T^ is seen to be 
under control for A^c = 4 and 6, as is the continuum Umit. For all Ac, the ratio Ae/Amax scales well. The numbers in brackets 
are the errors on the least significant digits, the first is statistical and the second systematic. 



related to the latent heat, we list the quantity Ae/Amax in the table. As one can see, this ratio has much better 
scaling properties, and the thermodynamic and continuum limits of the ratio are very well determined. 
The data collected for A't = 8 in Table H] is fitted extremely well by the form 

- 0.388(3) -HM. (9) 



where the numbers in brackets are the statistical errors on the last digit of the central values. Interestingly, the fit 
yields /S.e/ {(IaT^) = 0.014 ± 0.014 at Nc — 2, where there is a second order finite temperature transition, and hence 
Ae = 0. If one adds a term of 0{\/Nc), then the fit changes marginally: the limiting value for A'c — >■ oo is stable 
at the 3cr level. The coefficient of the ©(l/A^^) term changes by 16%, and the next correction term is marginal, its 
value being less than 10% of the total for Nc — 3. The extended series extrapolated to A'c = 2 is still consistent with 
vanishing latent heat of this theory. Although Nc — 2 must be the limit of the validity of the series expansion around 
1/Nc = 0, it seems to be well-behaved at iVc = 3. In agreement with this, a reliable value of the Nc ^ oo limit of 
Ae/ {cIaT'^) can be extracted. This is an example of successful strong scaling; Ae/ {(IaTc) is well described by just two 
terms of the series in 1/Nc even at small Nc. 

The series for the fourth root of the above quantity may be of interest, since v''Ae has mass dimension unity j3, Q- 
This gives 

(^)"*--™<')-^- 

The quality of the fit, as judged by the value of x^/DOF, is worse, but still within the limits of acceptability. This 
result above seems to have improved convergence properties around Nc — oo. However, on extending the fit to include 
the 0{1/Nc) term, we find that the correction terms are unstable against changes. The coefficient of the second term 
reduces to half its value, and the coefficient of the third term is 6-7 times larger. The values of these terms are 
comparable to each other for Nc = 3, opening the possibility that even higher order terms, or a resummation of 
the whole series, need to be taken into account. From the previous analysis it seems that the series expansion for 
Ae / {(IaTc) comes close to performing this resummation. We shall show later that other mass dimension four quantities 
such as p, e and s also have good strong scaling properties. 
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FIG. 4: Systematics of A/T'* in SU(4) and SU(6) lattice gauge theories. Both finite size and finite lattice spacing effects are 
small, except, possibly, very close to Tc. 
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FIG. 5: Testing a different scaling '2l'| of the conformal symmetry breaking measure: A/{dAT^T^), where dA = - 1, for 
Nc = 3, 4, and 6. 
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IV. CONFORMAL SYMMETRY BREAKING 



As discussed in section HIl A/T"* is easily calculated on the lattice. We have seen in the case of the latent heat, 
however, finite volume and cutoff effects need to be controlled. In the SU(3) gauge theory, it is known that A/T"' 
rises rapidly near and peaks at about LIT^ [l^. The rapid rise is, of course, dictated by the existence of a latent 
heat, but the shift in the peak away from Tc is not yet understood. We examine the Nc dependence of this peak. 

In Figure m we show the sensitivity of A/T"* to C and a. Except in the immediate vicinity of Tc there seems to 
be little sensitivity to C- The cutoff dependence is also insignificant, except in the vicinity of the peak, A,„ax. The 
good agreement between data obtained with a = 1/(6T) and 1/(8T) show that the results of the measurement with 
either lattice spacing can be taken to be an estimate of the continuum limit. We choose the conservative alternative 
of using a = 1/(8T) as a determination of the continuum results. The peak of A/T^ for SU(4) gauge theory is in the 
range 1.04 < T/Tc < 1.08. 

The systematics of A/T^ for SU(6) gauge theories is also shown in the same figure. The trends are very similar to 
those in SU(4). Results for different C agree very well. The approach to the continuum limit is also very similar to 
that discussed for SU(4). Again, in this case, we can take the results obtained with a = 1/{8T) to be an estimate of 
the continuum limit. In going from = 3 to 6, the peak in A/T'* moves closer to T^. 

It is of phenomenological interest to note that A/T'^ is not small even at 2Tc. In fact, as one can see in Figure SI 
one has 




2Te 



For Nc = 3 this implies A'^/^ ~ T = 2Tc. This is a natural scale, and therefore the theory is far from conformal. 

We end this section with the investigation of an intriguing observation made in j^^, [2l| : in the temperature range 
l.lTc < T < 4Tc, for the SU(3) pure glue plasma, A/T^ seems to be roughly constant. Phenomenological models 
of the gluon plasma have introduced mass scales and obtained such a behaviour [2^. We investigated this modified 
scaling behaviour at larger Nc (see Figure[5]). The dimensionless quantity A/ {(IaT'^Tc) for SU(3) gauge theory is seen 
to have little temperature dependence from just above Tc to about 4Tc. Unfortunately, the data for SU(3) theory is 
noisy at larger T. For SU(6) the error bars are smaller and one can observe that this quantity falls with T. This 
implies that the temperature dependence of A could be slower than . It would be interesting in future to expand 
the range of Nc and T in order to study this further. 



V. OTHER BULK THERMODYNAMIC QUANTITIES 

The pressure is calculated using the method outlined in ([6]). The integration requires interpolation of the measured 
points, and there could be a systematic error arising from this. We have estimated this error by comparing linear and 
quadratic interpolations. We found that point by point the error is small. Since the integration errors increase over 
the range of /? and reach a maximum at the highest /? that we use, it is sufficient to report the magnitude of that 
error compared to the statistical uncertainty. We can estimate the significance of the error by the i-statistic — 

t^\Il=I^= (12) 

where I2 is the integral estimates through a quadratic interpolation, Ii is the estimate using a linear interpolation and 
c^(/i,2) is the statistical error in the estimate of /i,2. This measure for SU(4) is 1.03 and for SU(6) it is 0.007. This 
source of error is therefore almost negligible. As a result, the systematic error is almost entirely due to the neglect of 
POj the pressure at the lowest temperature where the integration is started. 

Figure |6] shows the cutoff and volume dependence in the calculation of pressure for the SU(4) and SU(6) theories. 
The results are normalized by the known (asymptotic) finite cutoff correction for an ideal gluon gas, Psb, which 
was described earlier. Hence the pressure, so normalized, should go to unity. We find that finite volume effects are 
negligible. Finite lattice spacing effects also turn out to be negligible once we normalize the pressure by psB j3^ . 
These results indicate that it is safe to identify the continuum limit with the Nt = 8 measurements. Since the volume 
dependence is negligible, at each T we use the largest volume on which reliable results are available as an indication 
of the thermodynamic limit. 

At asymptotically large T the SU(A^c) gauge theory should go over to the ideal gluon gas. However, even at the 
highest temperatures which we have probed, i.e., T ~ 4Tc, the ratio p/psB is far from unity (see Figure |6l). Given 
p/T^ and A/T^ we obtain also the other bulk thermodynamic quantities: e/T^ and s/T^. These results are collected 
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FIG. 6: Study of systematics in the calculation of pressure for SU(4) and SU(6) gauge theory. Both finite size and finite lattice 
spacing effects seem to be under control. 
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FIG. 7: Results for thermodynamic quantities in SU(4) (left) and SU(6) (right) gauge theories. Shown are the pressure, energy 
density and entropy density, in units of pss. The ideal gas limits for these quantities should be 1, 3 and 4 respectively. 



in Figure [71 One sees clear deviations from the ideal gas limit in these two quantities as well. All three quantities 
also show a very slow rise throughout the measured range of T. Both e and s show a rapid jump near Tc, stronger 
for SU(6) than for SU(4). 
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FIG. 8: The first panel shows the energy density and pressure in SU(A'^c) pure gauge theories, normalized to the corresponding 
ideal gas values. The second panel shows A/T^, normalized by cIa = A^'^ — 1, for the same theories. 
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FIG. 9: The value of the 't Hooft coupling at the finite temperature transition, Ac, for different N^. The renormalized gauge 
coupling is evaluated at scale of 2ttTc, so that the result is independent of the lattice spacing. The boxes denote results obtained 
using the non-perturbative beta-function and the circles using the two-loop beta-function in the V-scheme. 



VI. Nc SCALING OF THERMODYNAMIC OBSERVABLES 



As discussed earlier, we distinguish between strong Nc scaling and 't Hooft scaling. The first is scaling with Nc 
of thermodynamic quantities at fixed T/Tc, and the second, the scaling with Nc at fixed A. We examine scaling by 
combining our results for the continuum limit of bulk thermodynamic quantities for Nc — 4 and 6, with a reanalysis 
of the older Nc = 3 data of [l3| using our techniques. 

Strong Nc scaling has been observed on the lattice in many contexts 0, S [IS [HI ■ The continuum limit of bulk 
thermodynamic quantities that we have extracted are also consistent with this limiting procedure, except for Nc = 3 
near Tc- In Figure [5] we show the energy density and pressure, each normalized by its ideal gas value, for Nc = 3, 4 
and 6. As noted earlier, there are clear deviations from the ideal gas behaviour, but the scaled quantities are almost 
independent of Nc- Since the ideal gas values scale with d^, one expects that at large T all bulk thermodynamic 
quantities should scale accordingly. 

We see that p/pss is independent of Nc over the whole temperature range above Tc within the accuracy of our 
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FIG. 10: The scaled entropy density, s/ssb, as a function of T/Tc and the 't Hooft coupling A — NcQb,- The gauge coupling is 
evaluated at the scale 2tiT using the non-perturbative beta-function for the appropriate theory. 



measurement. Such a statement is also true for the energy density except when T is close to Tc- Close to Tc the 
energy density does not scale as the ideal g 3/8, i.e., cLS 

dA, between Nc = i and larger values of Nc- This is, of course, 
a consequence of the fact that the latent heat density does not scale as dA for A^^ = 3 (see Table Consequently, 
A/r** also fails to scale with cLa in the vicinity of Tc- In Figure[8]we also show that the peak of A/ (T^d^i) is smaller 
for A^c = 3 than for A^c = 4 and 6 (the latter two are almost identical in value). We also see that the peak is rounded 
and shifted away from T^, at A'^ = 3. As discussed in Section Hill this could be due to finite-size effects, since N^ — i 
has a weaker first-order phase transition. With this exception, strong Nc scaling seems to work very well for bulk 
thermodynamic quantities; the sub-leading corrections are too small to be seen over the statistical errors. 

Next we turn to evidence for 't Hooft scaling. Since the critical point is known very precisely for several Q , we 
test this scaling using Ac = N cg\{'2,iTTc) [s^l- The change of Ac with A'c is very much larger than the statistical errors, 
and is seen using both the non-perturbative and the two-loop beta-functions, although it is somewhat larger with the 
former. The change is non-monotonic when the two-loop beta-function is used. In this context, we recall the result 
shown in the appendix: that the non-perturbative beta- function is preferred near Tc- Restricting ourselves to using 
only this leaves the three smallest values of A^c- 

Even with these three, it is clear from Figure [5] that the data do not fall on a straight line and hence a single 
correction term does not suffice. Using a second correction term, a description of the confinement-deconfinement 
phase boundary, i.e., the variation of Ac with is — 
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(13) 



Since the fit is linear in the parameters, the formal solution for minimization can be written down along with 
parameter errors. At A^c = 3 the values of the ©(l/A'f ) and ©(l/A^^) terms are comparable. Although they are small 
corrections to the leading, ©(A"^) term, this behaviour of the series could indicate that A"c = 3 lies near the radius of 
convergence of the series around A"c = oo . 

If this is so, then summing three terms of the series is not numerically accurate. However, with three pieces of data 
fitting a large number of terms is an ill-conditioned problem. As a result, one would do better to fit a resummation of 
the series, provided such a resummation has a small number of parameters. Unfortunately, there is no theory for the 
shape of the phase boundary. In its absence we try the usual trick of estimating a Fade resummation of the series. 
With three terms the best that we can try to do is to fit the lowest order expansion 



XciNc) = Ac(oo) + 



a/N^ 



1 + N^/N^ ■ 



(14) 



The best fit gives A^* ~ 4, roughly consistent with the series analysis. As a result, the fitted value of Ac(oo) is sensitive 
to the form of the remaining function, and cannot be reliably extracted using data for A^c near 4. It would be useful 
to improve the computations with Nc > 6 in order to extract this quantity with better accuracy. 

Next, we extend such a test to a bulk thermodynamic quantity. Our results for s/ssb as functions of T/Tc and 
A(27rT) are shown in Figure [TUl As we show in the figure, strong scaling holds with good precision since s/ssb is 
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FIG. 11: Approach to conformality in SU(A^c) gauge theories. The diagonal hne contains results from all possible conformal 
theories; the upper right end is the special case of a free massless theory. Also shown are the weak coupling results from [2^ 
(with increasing Nc the curves move up monotonically) . Since the lattice data lie closer to the weak-coupling results than 
to the diagonal line, there seems to be no window where a strongly coupled conformal theory describes SU(A'^c) pure gauge 
thermodynamics . 



almost independent of Nc down to the smallest temperatures that we have studied. For scaling at fixed A, convergence 
of the series is clearly bad for A ~ 9 or larger, because the theories for different Nc begin to drop out of the plasma 
phase; figure |9] shows that theories with smaller Nc drop to the confined phase at smaller A. At any fixed A in this 
region one has to go to Nc large enough that the theories are all in the same phase before one can observe good 
scaling with Nc- We also find that the convergence of the series in 1/iV^ is acceptable when A < 8.6 [33|. In the 
range 8.6 < A < 9, the physically interesting theory with A^'c = 3 is close to the radius of convergence of the series 
expansion, and the effect of the correction terms is large. 

Figure [10] also displays a comparison of the scaled entropy with predictions in a M — 4 supersymmetric Yang- Mills 
theory, computed |24j using the AdS/CFT correspondence — 

s 3 45 / IN 



Although one does not expect this computation to be valid in the realistic non-supersymmetric theories under inves- 
tigation, it is sometimes said to agree with lattice results. Here we show that the agreement is poor, except at the 
highest possible temperatures. It has been argued jTlj that more realistic AdS-QCD models should be used for such 
a comparison; the analysis of such models lies beyond the scope of this paper. 



VII. SUMMARY AND DISCUSSION 



In this paper we have studied the thermodynamics of the gluoA^c plasma, by numerical simulations of SU(4) and 
SU(6) gauge theories, and comparing them with a reanalysis of existing data [ij] for SU(3). Our focus is on taking 
the continuum and thermodynamic limits, by using multiple spatial volumes at each cutoff and by using significantly 
smaller cutoffs than used previously for the EOS of Nc > 3. As discussed in the introduction, in terms of statistics, 
lattice spacing and spatial volumes, this work brings the study of the thermodynamics of pure gauge SU(Afc) theories 
with small A'c > 3 at par with the state of the art for Nc — 3, while also throwing some new light on the SU(3) theory. 

We had shown in an earlier study of the deconfinement transition Q that observations made with finite lattice 
spacing could be continued to the continuum using the renormalization group equations. When measurements are 
made with lattice spacing a < l/(8Tc) we had found that the two- loop beta- function suffices. When the lattice spacing 
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0=1/ (6Tc) is used a nor-perturbative beta- function was introduced which could be used to continue the lattice results 
to the continuum limit 19|. In this paper we used these earlier results to obtain the continuum thermodynamics of 
pure gauge SU(4) and SU(6) theories. We also made a reanalysis of the older SU(3) data using this technique (as 
detailed in the appendix). 

One of our important results (see Section IIIip is an extraction of the latent heat of the deconfinement transition, 
Ae/T^'', for iV^ = 3, 4 and 6. We found that /\e/{dAT^), where dA = N'^ - 1, increases between A^c = 3 and larger Nc, 
indicating that the first-order transition grows stronger with increasing Nc- Our results are compatible with [8]. One 
expects stronger transitions to have smaller finite volume effects, and our observations support this notion. We found 
some scale breaking in the measurement of Ae/T^, and observed that Ae/Amax shows better scaling properties. We 
also saw that in the large Nc limit one has Ae/((i^T^) = 0.388 ± 0.003 (see eq. [S]). 

Further study of bulk thermodynamic quantities started with measurements of A/T'* (details are given in Section 
HV)) . We found good scaling of this quantity, and a reliable continuum limit, for T > Tc- For SU(3) some scale 
breaking is observed very close to Tc where the peak of this quantity lies. The cause remains obscure, although there 
are some indications that lead us to conjecture that this could be due to finite volume effects. Future studies are 
planned to understand this remaining ambiguity. In a range of temperature up to 4Tc we found that A cx [20l . [2ll | , 
or, possibly, slower. 

The pressure, p/T'*, was obtained using the so-called integral method (see Section |V|. Finite volume and lattice 
spacing effects in this measurement are under good control. We extracted the energy density, e/T'*, and the entropy 
density, s/T^, using these two primary measurements. In the whole range of T all these quantities lie substantially 
below the ideal gas values (see Figures [7] and [5]). Nevertheless, p/T'^, considered as a function of T/Tc, scales very well 
with dA [Sllll- Sub-leading corrections in Nc are hard to see at the level of accuracy we have reached (see Figure 
[5]). Similar scaling with dA is also seen for e/T"* and s/T^, except, possibly, in a small region near Tc- 

In Figu re [111 we present a plot of the normalized energy density, e/csB against the normalized pressure, p/psB 
following [2^. The diagonal line is the line of all conformal theories, with the ideal gluon gas being one special point 
on it. Weak coupling results [l^ lie below this. One sees that the lattice data for A''c = 3, 4 and 6 lie further below. 
Indeed the topology of these relations is such that the weak coupling predictions are always a better approximation 
to the lattice data than conformal theories. This figure gives clear evidence that there is no window of temperature 
in which the EOS can be described by a strongly-coupled conformal theory better than by weak coupling theory. 

Yet another reason for making high-precision measurements of bulk thermodynamics for Nc 7^ 3 is to understand 
the usefulness of the large- iVc limit (see Section |VI|. In agreement with previous results, we find that the strong 
scaling limit obtained by taking fixed T and Nc 00 works very well. Corrections in powers of 1/Nc are small, as a 
result of which the large- A^c results for the entropy density, for example, can be directly applied for Nc = 3 with about 
1% error. However, when the same data is analyzed as a function of the 't Hooft coupling, the finite Nc corrections 
are large for A > 8.6. The phase boundary for the large Nc theory, expanded in powers of i/Nc, seems to have a 
radius of convergence smaller than 1 /9 (see eq. [Til) . Since this observation has ramifications for all models of thermal 
QCD which proceed from the large- A^c approximation, including string-based models, we plan further measurements 
in the near future to explore the applicability of the 't Hooft limit. 

We thank Mikko Laine for providing us with the weak coupling results of Fig. 1111 and Rob Pisarski for comments. 
The computations were carried out on the workstation farm of the department of theoretical physics, TIFR and the 
Cray XI of the ILGTI. We thank Ajay Salve for technical support. 



Appendix A: The beta-function 



The lattice theory is cut off at a length scale of a. When a is small various quantities have a perturbation expansion 
in gn, the renormalized coupling determined at a scale na (where k can depend on the scheme). Then the derivative 



in eq. ([S]) can be written as 

dp dp dgn 



9 In a dgn, 9 In a ' 



(Al) 



where the last factor is the negative of the beta-function. Due to the fact that SU(iVc) gauge theories are asymptotically 
free, one expects that a weak-coupling determination of the beta-function should suffice when a is small enough. In 
it was shown that the two-loop beta-function is a sufficient description of the ffow of the renormalized coupling 
at the scale of a < l/(8Tc). Confidence in the efficacy of two-loop scaling is enhanced by the fact that the scheme 
dependence in the extraction of the QCD scale was small. 

Although the two- loop beta- function was insufficient to describe the ffow of the coupling at larger a, it was shown 
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FIG. 12: The step scaling function in SU(3) (left) and SU(4) (rig ht) gauge theory. The data points are for the location of the 
deconfinement transition. For SU(3) they are taken from [1^, [27| and for SU(4) from M. The non-perturbative beta-function 
of is tuned for use with a < ac(4). We have used a non-perturbative beta-function [a, [l^] which can be used for a < ac(6). 
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FIG. 13: A/T* for SU(3) gauge theory on 32^ x 8 lattices (left) and SU(4) on 30^ x 8 lattices (right). Shown are the analyses 
using various different beta- functions (based on the V-scheme). For SU(3), where the data is taken from [l^] we also show the 
analysis using the non-perturbative beta-function of that paper. 



that for a ~ ac(6), where ac{Nt) = l/{NtTc), a simple correction of the form 

with i?(a;) =e-^/Vi/2/3o^ (A2) 

suffices, where the two- loop beta- function is —fiog^ — Pig^ (for an alternative approach see [11]). The values of 
C2 were presented in Table (iv ) of In this paper the integration of the beta- function is started from the scale 
ac(8). In the calculation of section ITVl since we have taken data at lattice spacings as low as ac(6), we have used 
such a non-perturbatively corrected beta- function in the V-scheme (this differs from the conventions of 0]). Any 
non-perturbative beta-function will include finite lattice spacing corrections (s^l into the scaling, just as the above 
function does. Such corrections are non universal. 

A simplified version of the tests of scaling in [1, can be presented using the step-scaling function A/?(/3). This 
is the change in the bare coupling, A/3, required to reproduce the physics observed at a bare coupling /?, when the 
lattice spacing is doubled. If (3 is chosen to be the lattice coupling where the deconfinement transition is observed 
for a given Nt, then /? — A/? is the lattice coupling at the deconfinement transition when Nt is changed to Nt/ 2. In 
Figure [12] we show the result of using the beta- functions given above and the step scaling function given in [l3| . 

We also examined the sensitivity of the equation of state to the choice of the beta-function. Figure [13] displays 
resuhs for A/T^ obtained in SU(3) and SU(4) gauge theory, for lattices with a = 1/(8T). For the SU(4) theory, 
differences between the E-scheme and the V-scheme are statistically insignificant at all temperatures. Similarly, the 
difference between these and the non-perturbative beta-function of (|A2[) are also insignificant at these temperatures. 



aciNt)A - R 



\l3ogl{ac) 



16 



The results are similar for the SU(3) gauge theory, where we have re-analyzed the data of [l3|. The one-loop and 
the two- loop beta- functions give coincident results for T > LST^. In f28| a non-perturbative beta- function of the form 
in (|A2|) ) was used to describe scaling of Tc. The results of using this for A/T"' are shown. These three are close to 
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